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ABSTRACT 

Detection of signals from a possible extrasolar technological civilization is one of the challenging 
efforts of science. In this work, we propose using natural telescopes made of single or binary 
gravitational leasing systems to magnify leakage of electromagnetic signals from a remote planet 
harbours an Extra Terrestrial Intelligent (ETI) technology. The gravitational microlensing surveys 
are monitoring a large area of Galactic bulge for searching microlensing events and they find more 
than 2000 events per year. These lenses are capable of playing the role of natural telescopes and 
in some occasions they can magnify radio band signals from the planets orbiting around the source 
stars in gravitational microlensing systems. Assuming that frequency of electromagnetic waves used 
for telecommunication in ETIs is similar to ours, we propose follow-up observation of microlensing 
events with radio telescopes such as Square Kilometre Array (SKA), Low Frequency Demonstrators 
(LED) and Mileura Wide-Field Array (MWA). Amplifying signals from the leakage of broadcasting 
by an Earth-like civilizations will allow us to detect them up to center of Milky Way galaxy. Our 
analysis shows that in binary microlensing systems, the probability of amplification of signals from 
ETIs is more than that in single microlensing events. Finally we propose target of opportunity mode 
for follow-up observations of binary microlensing events with SKA as a new observational program 
for searching ETIs. Using the optimistic values for the factors of Drake equation provides detection 
of about one event per year. 


1. INTRODUCTION 

Are we alone in the Universe ? This is one of the deepest questions of mankind and in recent years there have been a 
lot of efforts to answer this question. The Search for Extra Terrestrial Intelligence (SETI) is one of the projects trying 
to answer this question and since 1960s no confirmed signal from extraterrestrial life has been detected (Wilson 2001). 
Another important effort started with direct and indirect observations of extra-solar planets after the discovery of first 
exoplanet by Mayor & Queloz (1995). Nowadays, projects like Kepler increased the number of extrasolar planets to 
more than one thousand candidates (Batalha el ai. 2014) and an important question still is, are there any intelligent 
life exists in one of the discovered planets ? 

In order to detect Extra Terrestrial Intelligent (ETI), we need to listen leakages from the telecommunication signals 
of advanced technology civilizations. The transmitting signals can be either particles or photons. The properties of 
transmitting particle has been discussed extensively in Oliver & Billingham (1996) with the following properties of (a) 
having minimum energy per quantum for transmitting particle, (b) large velocity for particle (c) particles easily generate 
and focused and finally (d) particles have not easily absorbed by the interstellar medium. The best candidate that 
satisfies these four properties is the electromagnetic waves in radio or microwave wavelengths . The main challenge of 
detection of signals from an advanced civilization is that in the large astronomical scales signals are faded out. Also, in 
the radio waves for frequencies smaller than lOOMHz (larger than 3 meter), the free-free absorption become important 
and for frequencies larger than ITHz (smaller than 300 micro-meter) absorption becomes increasingly important. 

Assuming that the electromagnetic signals are transmitted by an advance civilization, these signals should have the 
specific property of violating the natural transmission mechanism to be distinguished from the natural sources (Tarter 
2001). Moreover, we expect that transmission by the intelligent life has to be done in narrow band by the artificial 
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Service 

(1) 

Freq. 

(2) (MHz) 

Transmitters 

(3) (No.) 

Max. Power 

(4) perTr. (W) 

Bandwidth 

(5) (Hz) 

Power 

(6) (W) 

Power/Hz 
(7) (W/Hz) 

Military 

~400 

10 

2 X 10“ 

10® 

2 X 10“ 

2 X 10® 

FM 

88-108 

9000 

4 X 10® 

0.1 

4 X lO’’ 

4 X 10® 

TV 

40-850 

2000 

5 X 10® 

0.1 

10® 

10“ 


Table 1. Characteristics of the strongest radio transmitters on the Earth, adopted from Sullivan et al. (1978) and 
Loeb & Zaldarriaga (2007) .The power of military transmitter is corrected based on Loeb & Zaldarriaga (2007). The first 
column is the type of service, second column is the frequency range of transmission, third column is the number of transmitters, 
forth column is the maximum power per transmitter, fifth column is the bandwidth, sixth column is the total power of transmit¬ 
ters and the seventh column is the transmitter power per Hertz. Moreover, Arecibo telescope has four radar that transmitting 
with the power of 2 x lO’^^ W (Tarter 2001). 


sources with the minimum amount of AA/A, similar to what our civilization is doing for transmitting signals. With our 
primitive technology, we can produce AA/A ^ 10“^^ compare to nature that is around 10“'^ to 10“^ in astrophysical 
systems and water masers with limiting value of 10“®. Assuming that Earth-like advanced civilizations using the same 
range of electromagnetic waves for telecommunications that we use, the present radio telescopes can eavesdrop their 
radio transmissions (Loeb & Zaldarriaga 2007). 

The aim of this work is using gravitational microlensing as natural amplifiers to magnify leakage of signals from 
remote planets at the Galactic scales. The gravitational lenses can magnify electromagnetic waves from distant sources 
by converging light beams from the source star. This light magnification can be done either by single or binary lenses. 
At the present time, there are two main observational campaigns of OGLE and MOA that use wide field telecopies to 
monitor centre of Galaxy for gravitational microlensing detections. In 2014, OGLE could detect over 2000 microlensing 
events (Sumi et al. 2014; Gaudi 2012). The next generation of microlensing surveys as Korean Microlensing Telescope 
Network (KMTNet) is going to cover a larger area and discover more microlensing events than the present experiments 
(Park et al. 2012). 

The main propose of this work is the follow-up observation of microlensing events with radio telescopes to detect 
transmitting signals from planets hosting ETI and orbiting around the source stars in the microlensing systems. During 
the gravitational microlensing, electromagnetic waves independent of their wave-lengths are magnified from both the 
parent source star and a planet orbiting around it. However, since the source star and the associated planet are 
separated with the distance in the order of astronomical unit, their magnifications do not peaks at the same time. 
This method increases the depth of SETI observation from few parsecs up to the centre of Galaxy. This technique 
in the visible band also has been proposed for detecting hot Jupiters by illuminating planets around the parent stars 
(Sajadian & Rahvar 2010). 

The organization of this paper is as follows: In section (2), we introduce the radio sources on an Earth-like planet 
similar to our transmitters on the Earth and study the possibility of detecting these sources with the present instruments 
that are using in the cosmological observations. In section (.->), we introduce the basic formalism of gravitational 
microlensing as a natural amplifier of signals from an Earth-like planet by single and binary leasing systems. In 
Section ( ), we perform a Monte-Carlo simulation and estimate the number of high magnification microlensing events 
of planets by follow-up radio observations for both single and binary leasing systems. In this section we also investigate 
the wave optics effects of gravitational microlensing in our observation and the suppression of very high magnified 
singles in the radio light curve. Details about this effect is given in appendix (A). In Section (.4) after introducing 
Drake equation, we investigate the number of expected candidates with signals from the planets with intelligent life. 
The conclusion is given in section (6). 

2. RADIO SOURCES FROM AN EARTH-LIKE CIVILIZATION 

Let us imaging an Earth-Like civilization that is producing at the same frequency rang and the same power of 
electromagnetic waves that we are producing on the Earth. We note that according to the classification of advanced 
civilizations by Kardashev, these planets may produce stronger transmissions compare to the Earth (Kardashev 1964) 
or due to improving fibre optics communication on their planet, they produce less man-made noises compare to ours. 
The transmitted signals might be either in beacon mode or unintentional leakage of signals. Here we use the later 
mode of transmission in our study (Loeb & Zaldarriaga 2007). As we discussed before, the artificial signals are narrow 
band and encoded signals contains more information compare to the natural signals. The complexity of signals from 
intelligent civilization can be measured by various mathematical methods. 
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Figure 1. The logarithm of Power per area (in TT/m^) from a Military-like source adapted from Table (1) transmitted from 
an Earth-like intelligent life as a function of distance from the source (in parsec). The horizontal dashed line (i.e. P/S = 
/m^) is the limiting value that SKA can observe signals, means that the source is detectable up to distance of D ~ 1.26 

kpc. 


Different types of radio sources as a target for intelligent civilization has been discussed in Sullivan et al. (1978). Ac¬ 
cording to Table (1) which is adapted from Sullivan et al. (1978) and Loeb & Zaldarriaga (2007), military radars with 
maximum power per transmitter {Pmax — 2 x 10^^ W) is the strongest sources that radiate radio-waves isotropically. 
Television and FM stations are the second and third ranks in the power of radiation. Moreover, Arecibo Observatory 
uses transmitters that use them occasionally for ETI communications. This observatory has four radar transmitters, 
with effective isotropic radiated powers of 20 TW at 2380 MHz, 2.5 TW (pulse peak) at 430 MHz, 300 MW at 47 
MHz, and 6 MW at 8 MHz (Tarter 2001). Assuming an Earth-like planet in our Galaxy, SETI project can investigate 
the following characteristics to search for the signals of intelligent life, (a) short bandwidth transmission, (b) frequency 
and intensity modulation due to rotation of planet around its axis and around the parent star and (c) measurement 
of complexity in data encoded in the signals. 

In recent years the cosmological experiments are designed to observe 21cm emissions of HII regions from redshift 
of z ~ 6 — 12, corresponds to 200 — 400 MHz (Becker et al. 2001). These cosmological observations are aiming to 
observe the reionization history of Universe. Fortunately, this range of frequency overlaps with the frequency range of 
telecommunications on the Earth (see Table 1), and it can be used for observation of Earth-like planets in our Galaxy 
(if intelligent life uses the same range of frequency). While military transmitters are the strongest radio sources, TV 
and FM radio sources with 10 — 100 kW can outshine Sun up to almost 100 times, means that radio emission from 
parent star does not blend with signals from an Earth-like planet. So the main limitation on the observation is the 
threshold of signal to noise ratio of receivers that we are using on the Earth. The signal to noise increases by increasing 
the square root of exposure time of observation, however exposures longer than the spinning period of planet or orbital 
period of planet around the parent star causes broadening of emission signals and erasing the modulation effect. On 
the other hand for long exposures, we can not detect complex signals from the radio flux by accumulating signals. 
Here we study the sensitivity of MWA-LFD and SKA for detection of signals from the intelligent life. 

For the case of MWA-LFD, the point source sensitivity (PPS) at 200MHz is ^ 


Fpss = 0.4 mjy 


Av fobs \ 
8kHz month J 


( 1 ) 


where tabs is the total exposure time and Au is the bandwidth. For a source with a minimum power of Pmin located 


^ The value of FPSS for other frequencies can be found at http://web.haystack.mit.edu/arrays/MWA/LFD 
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at the distance of D and emitting with the bandwidth of A^, the threshold of observability from equation (1) is 


= 1.99 X X 


D 


100 pc 


1 kHz ) ^ 


^obs 

10 min 


- 1/2 


(2) 


Let us assume that we perform observations with the exposure time of 10 min. This strategy of observation enable us 
to monitor modulation of frequency within few hours from the spin of an Earth like planet. Comparing the minimum 
flux of detection with the power of military transmitters in Table (1) as well as Arecibo transmitter, we can detect 
them up to lOpc and within this volume there are almost thousand stars available. This volume can be increased with 
increasing the exposure time of observation, however we will miss the sensitivity to the modulation. 

The other important telescope in our desired range of frequency is the Square Kilo-meter Array (SKA) telescope. This 
telescope also is designed for astrophysical and cosmological studies. It has the threshold sensitivity of 10“^® Wjrr? 
and is able to monitor sources with high cadence. This will enable us to monitor time variation of source flux and 
search for possible meaningful signals from an Earth-like planet. By using narrow bandwidth of detection, it is possible 
to improve the threshold sensitivity to 10“^® Wjw? (Lazio et al. 2004). Using the strongest transmitter on an Earth¬ 
like planet from Table 1, we can put limit on distance that can be detected with SKA. The result from Figure (1) is 
Dq ~ 1.26 kpc. 

In the following sections we will introduce gravitational microlensing as the natural amplifiers of signals from remote 
planets. We will provide an observational strategy based on observational mode of microlensing surveys and estimate 
the number of planets that are illuminated with the gravitational microlensing in radio band. 


3. GRAVITATIONAL LENSING AS NATURAL TELESGOPES 


The gravitational lensing of a star by another star in the Milky Way galaxy, for the first time is introduced by 
Einstein in (1936) and since angular resolution of an optical telescope is less than the Einstein ring, Einstein noted 
in his paper that ” There is no great ehance of observing this phenomenon^'’ Einstein (1936). On the other hand the 
probability of detecting gravitational lensing inside the Galaxy is of the order of 10“^. After five decades, Paczyhski in 
1986 revisited the question of gravitational microlensing and estimate the observability of this phenomenon inside the 
Galaxy and propose the application of it for searching compact dark matter objects inside the Galactic halo Paczynski 
(1986). 

In the gravitational lensing inside the Milky Way, while the angular separation between the images is small and by 
the ground-based telescopes images are unresolvable, the magnification of light from the source star during the transit 
of lens is a measurable quantity. The magnification depends on the angular separation between the source and lens 
stars. The overall magnification from single lensing is given by 


A{t) = 


2 -I- u{tY 
M(t)\/4 -I- u(t)2 ’ 


(3) 


where u is the impact parameter normalized to the Einstein ring. The impact parameter is a dynamical parameter 
and for a single lens changes by time as w(t)® = Uq + {t — tof' h’E where uq is the closest distance from the lens, to is 
the moment of maximum magnification and tE is the Einstein crossing time (Paczynski 1986; Rahvar 2015). 

We imagine that in addition to the source star, there is a planet with intelligent life orbiting around the source star. 
The projected orbit of planet on the lens plane and relative distance of these three objects is shown in Figure (2). 
We take r as the orbital radius of planet projected on the lens plane, normalized to the Einstein radius of lens (i.e. 
f = w as the angular velocity of planet and 4> as the initial phase of planet. The equation of motion of planet 

in this coordinate system is given by 


Xp = f cos{ujt -I- (/)), 

L), = f sin(a;<-I-<)>) sin,5, (4) 


where 90 — /3 is the angle between the Y-axis and the orbital plane and the orbital plane goes through the X-axis. The 
parent star is located at the centre of coordinate system and the trajectory of lens is given by 


Xl = { 

Yl = { 


t — tp s 
tE 

t — tp s 
tE 


cos a — uq sin a, 
sin a -I- Uq cos a. 


(5) 


where a is the angle between the lens trajectory and semi-major axis and uq is the minimum impact parameter 
compare to the source star. The distance between the projected position of planet on the lens plane and lens, upL, is 
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Figure 2. The projection of source star (represented by star) and the orbital path of planet (elliptical) on the lens plane. The 
lens (triangle) moves on a straight line on the lens plane. The minimum distance of lens to the source star is given by uo and f 
is the projected orbital radius of planet on the lens plane, normalized to the Einstein radius. 




Figure 3. Typical magnification of flux from a planet as a function of time (normalized to the Einstein crossing time). Left 
panel represents microlensing light curve by a single lensing and right panel by a binary lensing. 


calculated based on relative distance (projected on the lens plane) between the two objects as follows: 

n 2 


WpL(i; to, Uo, tE,a, u}, r, (3) = 


(^—cos a — Uo sina — fcos(ujt + </>) 
ts 


(-^) sin a-\- Uo cos a — f sin(a;t + (f)) sin /3 

ts 


( 6 ) 


We can calculate the magnification of planet’s light by substituting this impact parameter in magnification formula 
for single lensing in equation (.->). Here the impact parameter, upL(t), depends on various orbital parameters of planet 
as f, w and /3 and the lens parameters as tp, uo and a. 

Figure (••■) on the left panel, represents the microlensing magnification factor of a planet orbiting around the parent 
star by single lensing. We note that the observation of microlensing events in short and long wavelengths has a main 
difference where for the long wavelengths and a point-like source, we have to take into account the wave optics features 
(Mehrabi & Rahvar 2013). In this case images on the lens plane plays the role of slits as in the Young experiment 
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Figure 4. OGLE field for microlensing events in the direction of Galactic Bulge. This figure is adapted from Udalski et al. 
(2015). 

and result is producing diffraction pattern on the observer position. The observational consequence is reducing the 
magnification factor of leasing and avoiding mathematical singularities in the geometric optics. This effect is discussed 
in Appendix (A). 

In the case of binary leasing, the formalism of leasing is complicated than single leasing (Witt & Mao 1995). In this 
case binary lens produces caustic lines (or singular lines) on the source plane. Whenever source crosses these lines, the 
magnification gets infinite. However in reality, due to finite size of source star and interference of coherent light, the 
magnification is a finite value. Figure (o), on the right panel shows a typical light curve of a binary lens magnifying 
signals from a planet orbiting a source star. Here we consider the geometric optics in calculation of magnification 
(Doininik 2007). 

In the next section we estimate the number of planets that their signals from intelligent civilizations can be illumi¬ 
nated by single or binary leasing up to the observational threshold of SKA telescope. 

4. STATISTICS OF PLANET ILLUMINATION 

In this section we estimate the number of Earth like planets that their signals can be amplified with gravitational 
microlensing (as a natural telescope). The observation of microlensing events are done in two steps of (a) the survey 
mode, where survey groups are monitoring stars in the direction of Galactic Bulge and by ongoing analysis of light 
curves, alert the microlensing events, (b) The follow-up mode, where telescopes around the globe observe microlensing 
candidates with high cadence and better photometric precision. The map of survey for OGLE experiment is shown in 
Figure (4) towards the direction of Galactic Bulge . In this section, we will estimate the statistics of high amplification 
of Earth-like planets by single and binary leasing systems. 

4.1. Single Lens 

One of the crucial parameters in simulation of this effect in single and double lensing is determination of f as the 
projected orbit of planet on the lens plane, normalized to the Einstein radius as follows 

- ^ rg(Ms) ( Ml \ ^ 

^ 5.7 AU \0.5Mq) \8kpcJ VI-a; 

where x = Dl/Ds represents the ratio of observer-lens to the observer-source distances, Ml is the mass of lens and 
rn is the orbital radius of a planet assuming that it is located at the habitable zone of a source star. We note that 
rn is a function of mass and the age of source star that planet orbits around it. Eor the main sequence stars, it is 
given by rniMs) = (0.8 — 1.5) x [Mg/M qY AU (Kasting et al. 1993; Koppapapu ct al. 2013). For the red giants, 
while these type of stars are unstable, however for a solar type star stay at the first stages of its post main-sequence 
evolution, the temporal transit of the habitable zone is estimated to be of several 10® years at 2 AU and around 10® 
years at 9 AU (Lopez et al. 2005). 

Fortunately from studying the color-magnitude distribution of source stars in OGLE HI microlensing events, we 
can estimate the statistics of stars in the main-sequence and in the red giant area. For OGLE HI in the direction 
of Galactic bulge, we adapt statistics of microlensing events from the observations in the period of years 2001-2009. 
Analyzing the distribution of stars in Figure (5), the fraction of red-clump source stars to the overall number of stars is 
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Color—magnitude diagram for OGLE-III events 
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Figure 5. Color-magnitude of source stars of microlensing events in OGLE III survey from year 2001-2009 in the direction of 
Galactic bulge. The magnitude of source stars are unblended after imposing the blending factor from the best fit. The top-right 
area of this diagram represents the red-clump source stars that comprises 16% of all microlensing candidates. This figure is 
adapted from Wyrzykowski et al. (2015). 


16 percent. We use the same combination for the population of source stars in our Monte-Carlo simulation, associating 
a habitable zone to the source stars. 

Taking into account that the limiting distance for detection of a military transmitter from an Earth-like planet 
with the SKA telescope is Dq = 1.26 kpc (for details refer to Section 2), the minimum magnification factor from 
microlensing is related to detectability with SKA telescope at the distance of D as follows 


A 


min 


> 


Do) 


2 


( 8 ) 


For the case of microlensing events in the direction of Galactic Bulge, if we assume for instance a source star located 
at Ds — 7 kpc, the minimum magnification for the planet detection is Amin > 30. For the microlensing events in the 
direction of Galactic Bulge, the source stars are not located at a fixed distance, rather they are distributed along our 
line of sight from the position of the observer to the centre of Galaxy. The main sequence stars are mainly belong to 
the disk and red clumps belong to the Galactic bulge. 

In order to estimate the fraction of events with a magnification that satisfy the condition in expression (8), we perform 
a Monte-Carlo simulation where the geometrical parameters of lens is chosen uniformly in the range of uq G [0,1], 
a G [0, 27r], /3 G [0, tt] and 4> G [0, 2tt]. The other physical parameters as the Einstein crossing time tE and r are taken 
according to the distribution of matter, dynamics of Galaxy and the mass function of lens stars in the Galaxy. The 
mass function including main sequence, brown dwarfs and compact stars are given with ^(logm) = dn/dlogm, as 
follows Chabrier (2003) 


^(log to) = 0.093 X exp 
^(logm) = 0.041 X to“^'^® 


(log TO — logO, 2) 
2 X (0.55)2 

, TO > TO0 


21 


TO < TO0 


(9) 


On the other hand, for simulating the source of the microlensing events, we adapt a complete sample from the Hipparcos 
catalog Turon et al. (1995) and enhance the contribution of red clumps to be comparable with the observation in the 
Galactic centre. 

For the distribution of matter in the Galaxy, the density of lenses composed of density of disk and bulge is described 
with a thin disk and a central bar structure. The disk is exponential model adapted from Besancon model Robin et al. 
(2003) 


Pd{r,z) = po 




age < 0.15 Gyr, 


exp(—(0.5^ 


)i/2)_exp(-(0.52+ 






( 10 ) 


Pd{r,z) = po 


age > 0.15 Gyr, 


( 11 ) 
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where = 5.0 kpc and /i_r_ = 3.0 for age< 0.15 Gyr and = 2.53 kpc and = 1.32 for age> 0.15 Gyr, 
+ (z/e)^. e and po for different ages of stars are given in detail in Robin et al. (2003). 

The bar is described in a cartesian frame positioned at the Galactic center with the major axis x tilted by $ = 13 
degree with respect to the Galactic center-Sun direction. The number density of starts in the bulge also is given by 
Robin et al. (2003) 

nB{x,y,z)=noe~'^^^'^, < Rc ( 12 ) 

UB^x^y, z) = nQe~'^ / 2 g- 2 (-;/a; 2 +y 2 _/j^) ^ ^-y"^ > Rc (13) 


where = 


[(!)" + (!)’]’ + ({)' 


11/2 


no = 13.7pc a = 1.59 kpc, h 


0.424 kpc, c = 0.424 kpc and Rc = 2.54 


kpc are the scale length factors. 

For the dynamics of Galaxy we apply separately the dynamics of stars in the disk and bulge. The global rotation of 
the disk is given as a function of the galactocentric distance by 


Vrot{r) 


Vrot,Q X 1.00762 



0.0394 


0.00712 


(14) 


where Vrot,Q = 220km/s Brand & Blitz (1993). The peculiar velocity of the disk stars is described by an anisotropic 
Gaussian distribution with the following radial, tangential and perpendicular velocity dispersions of cr^ = 34 km/s, 
ag = 28 km/and cr^ = 20 km/s Rahal et al. (2009). The distribution of transverse speed of stars in the bulge also is 
given by 


frivx) = — -?;_Lexp 

^ bulge 



(15) 


where a bulge — 110 km/s. 

The relative transverse speed of source and lens on the lens plane is given by vj_ = a:(vs — V£;)j_ — (v^ — V£;)j_ where 
all the vectors define in two dimension parallel to the lens plane and perpendicluar to the line of sight. The rate of 
events, taking into account that each lens spans a tube with the diameter of 2Re and the length of w_l x Tabs is given 

by 

dr oc ^Mlx(\ - x)p^,{x)^(\ogML)v±fs{vs)fL{vL)dxdvsdvLd\ogML, (16) 


where p*(a;) is the stellar density of Galaxy, ^ is the mass function of stars in logarithmic scale, fs{vs) and JhivL) 
are the velocity distribution functions of source and lens in two dimension. The important observable of microlensing 
events is the duration of event which is given by Ie and detectability of an event directly is related to this parameter. 
In our simulation, we adapt the detection efficiency function (i.e. e(tE)) from the OGLE survey Wyrzykowski et al. 
(2015). 

For single-lens microlensing events in the direction of Galactic Bulge, assuming a planet orbiting at the habitable zone 
of source star, almost 2.5% of overall observed microlensing events satisfies the condition of expression (8). Figure (6) 
represents the distribution of microlensing events in term of distance of source stars (i.e. Ds) for overall microlensing 
(in solid line) and the distribution of microlensing events satisfy the condition of expression (8) (in dashed line). The 
average distance of observer to the source stars for all the events is 8.36 kpc and for those with positive signal from 
planet is 7.49 kpc. While the observation of microlensing events with nearby sources improve receiving strong signals 
from the planets, estimation of distance of source star in the ongoing microlensing events is a challenging problem. 

To estimate the distance of source star, the color-magnitude of stars around the source star within the radius of 60" 
is measured (as an example see Street et al. (2016)). Assuming that all the stars in the field are affected by the same 
amount of extinction, using the position of center of red clumps in the color-magnitude diagram from one hand and 
the position of red clump stars in the Hipparcos catalog (Paezyhski & Stanck 1998), we can calculate the extinction 
of the held. Then the dereddened color and magnitude of the source star (i.e. Ao, (fo — ^)so) with measuring its offset 
with the red clump can be calculated. The result would be an estimation from the distance modulus of the source star 
as well as the angular size of the source star. Due to uncertainty in the distance measurement, we propose performing 
the radio follow-up observation for all the microlensing candidates regardless of their source distances. 

In this part, we have introduced the gravitational microlensing as a natural ampliher to enhance signals from planets 
and enable us to detect them up to distances of 10 Kpc. In another word, the depth of our observation with SKA 
could be extended up to the center of Galaxy. With the present technology, microlensing surveys observe almost 2000 
microlensing event per year. Regarding that the efficiency for planet illumination from single lensing in the Monte-Garlo 
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Figure 6. The distribution of distance of source stars (i.e. Dg) for overall microlensing events from the Monte-Carlo simulation 
in solid line and the distribution of microlensing events satisfy the condition of (8) with 2.5% of overall events in the dashed 
line. The distribution is scaled in logarithmic scale. Magnifying signals from planets in the microlensing event, we can detect 
them up to the centre of Galaxy 


simulation is about 2.5%, we would expect to detect almost 50 high magnification of planets around the microlensing 
source stars. For these events the average value of orbital size of habitable planets normalized to the Einstein radius 
is < f >= 0.18, compare to < f >= 0.19 which is the average value of this parameter for all the microlensing events 
regardless of planet detection. In the next section we will estimate the number of high magnification planets by the 
binary lensing. The probability of these planes hosting intelligent communicating life will be discussed later. 

4.2. Binary lensing 

In the binary lensing, the lensing system composed of two objects. These objects can be either (a) a binary star 
or (b) a star with its companion planet or (c) any other combination of sufficiently high mass objects, such as black 
holes, brown dwarfs, etc. Increasing the number of lenses from one to two or more objects increases dramatically the 
complication of the lensing equation. One of the features of having multiple lensing is the caustic lines on the source 
plane. Once a source crosses these lines, in theory we have infinite magnification and in the observation, we can get 
sparks of light from the source object. In the case of binary lensing, either source star or planet orbiting around it 
or both objects can cross the caustic lines. The advantage of using binary lensing for detecting planet is that with 
an early warning system, from the deviation of light curve from the single lensing, we can identify an ongoing binary 
microlensing event before the light curve peaks. We will show that the probability of high magnification of planet in 
binary lensing is more than a typical single lens. 

Investigating the early warning database of OLGE III surveys during years 2002-2003 results in detection of 15 
binary microlensing events out of overall 851 microlensing candidates. In another word, 1.76% of microlensing events 
in this catalog are binary lenses. Jaroszynski ct al. (200-1) found the best fit to the light curve and could reconstruct 
the caustic lines on the source plane or the projection of these lines on the lens plane. 

In what follows we adapt the caustic structure of binary lensing from the best fit. The stellar type of source stars are 
estimated from the position of stars in the color-magnitude diagram and for each star we can calculate from equation 
(7) a habitable zone for the planet orbiting around the source star. Then we calculate the probability of crossing of 
a planet orbiting around the source star by the caustic lines of a binary lens. We adapt the distribution of distance 
between the lens stars as in the observation sample that is shown in Figure ( ). For comparing the size of orbital radius 
of planets in the habitable zone with the caustic lines and probability of caustic crossing, we repet the Monte-Carlo 
simulation as in the single lensing. 

Taking into account that source stars with positive planet signal is located almost at < Ds >= 7.49 kpc, this 
distance corresponds to the distance modulus of /r = 14.37. On the other hand from studying red clumps in the 
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Figure 7. Distribution of distance between lenses of 15 binary micolensing events from OGLE 2002-2003 microlensing list 
(Jaroszynski et al. 2004). The distance between the lenses are normalized to the Einstein radius of overall mass of lenses. 



Figure 8. Distribution of f for microlensing events in the direction of Galactic bulge from the Monte-Garlo simulation. The 
average value of f is 0.19. Y-axis is in the logarithmic scale. 

Galactic Bulge, the average extinction of stars is < Ai >= 1.90 (Paczyhski & Stanek 1998). Then the average overall 
distance modulus with taking into accout the extinction for the direction of Bulge events would be ^ = 16.27. We 
use Table 1 from laroszynski et al. (2004) and obtain the absolute magnitude of stars in 1-band ranges from —0.64 to 
3.00. Comparing with the Hipparcos absolute magnitude (Turon et al. 1995), we can conclude that these stars are in 
the bright part of stellar population. In the gravitational microlensing experiment, we have a bias for observing bright 
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stars as the source of microlensing events in direction of Galactic bulge. Using the habitable zone for the orbital radius 
of source stars from equation (7) , we plot in Figure (8) the distribution of r where the microlensing events have been 
triggered with the microlensing surveys. The average value of this parameter is obtained as < f >= 0.19. 

Now we compare the size of r with the distance between the binary lenses (i.e. d). The average value of d from 
the observations of 15 binary lenses in OGLE 2002-2003 season is < d >= 1.26. The comparison of these two length 
scales provides the probability of crossing of planet’s orbit with the caustic lines of binary lenses. We note that due to 
binary feature of microlensing light curve from the source star, source star crosses the caustic lines. Also for the case 
of equal mass binary lenses, in the intermediate regime where the distance of lenses is in the range of l/v^ < d < 2, 
the caustic lines are connected bar shape on the lens plane (Dominik 1999) where comparing the intermediate regime 
for the caustic structure with Figure ( ) means that almost all the cases of binary lenses are in this regime. 

In order to estimate the probability of caustic crossing of a planet, we take a simple faced-on orbital configuration 
of planet around the source star and consider the relative motion of the projection of the planet’s orbit on the lens 
plane with respect to the binary lens. Since r <C d, the probability of caustic crossing of planet can be estimated by 


P ~ 


d — f 
d 


(17) 


where using the numerical values for f and d results in P ~ 0.85. Taking into account that 1.75% of microlensing 
events in our sample are binary lenses, then the probability of caustic crossing of a habitable planet in the binary 
lenses would be 1.5% of all the microlensing events. 

All our calculation is based on geometric optics, however having a point mass source with the longer wavelengths 
emission from the source produces wave optics features from the microlensing magnification and result is suppressing 
the magnification and avoiding the singularities. In the appendix (A) detailed calculation on the wave optics and its 
effect on the result of this section is discussed. We show that even using the gravitational leasing in the wave optics 
formalism doesn’t change the results of this section. 


5. STATISTICS OF INTELLIGENT LIFE DETECTION IN FOLLOW-UP MICROLENSING OBSERVATIONS 

We have seen in the previous section that the number of planet illumination from the the single lensing is almost one 
order of magnitude more than that in the binary lensing. However, from the the analysis of Monte-Carlo simulation, 
for almost 85 percent of the binary lensing there will be caustic crossing of planet, means that there will be high 
magnification by illuminating signals from the planet. Using binary lensing events for the follow-up observations has 
two practical advantages of, (a) in the binary microlensing events, from the deviation of light curve, we can identify the 
binary microlensing systems in advance where the probability of illumination of habitable planet is almost 85 percent, 
(b) since the probability of caustic crossing of planet is high, SKA telescope can be used in target of opportunity mode 
for the follow-up observation. 

In what follows we estimate of number of intelligent communication life in our observation method. We use the 
Drake equation which provides the probability of existing a Communicating Intelligent Life (CIL) in our Galaxy. This 
equation is given by 

Ncil = NsfpnEfLfifcL/Ls, (18) 

where Ns is the number of stars in our galaxy, fp is the fraction of stars with planet, np number of Earth like planet 
in the habitable zone, fp is the fraction of planets that life is developed, // is the fraction of those planets that evolve 
intelligent life and fc is the fraction of those intelligent life that can have radio communication. L is the time scale 
that intelligent civilization exists and Ls is the overall age of star in the stable position. 

There are optimistic and pessimistic estimations for different factors in Drake equation. We adapt updated estima¬ 
tions for each element. Observations of extrasolar planets with microlensing reveals that stars can have more than 
1.6 Earth like planet means that fp is at least one (Cassan et al. 2012). On the other hand recent analyses of the 
Kepler satellite showed that about 20% of all Sun-like stars have Earth-sized planets orbiting within the habitable 
zone (Petigura ct al. 2013) (i.e. np = 0.2). Also HARPS team estimated that almost 30% of Solar type stars have 
planet with an Earth mass or Super earth mass (Pepc ct al. 2011). If we assume that in these planets life eventually 
develop, means that the probability of formation of life is one, then we can write fpnpfp ~ 0.2. 

The later factors of fi, fc and L is unclear since only we have data from one case that exists on the Earth. We 
don’t know what fraction of life supporting planets eventually develop to an intelligent life. Does it a fixed point of 
life evolution and with the probability of unity we will have an intelligent life or the probability of intelligent life is 
rare. Does emerging intelligent life need a long time in the orders of billions of years or it can happen in a short time 
scale? Also how long an intelligent life can exist. Would it disappear by self destruction or it may long live with 
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spreading over other planets, even in the form of artificial intelligent life. We recall the multiplications of these factors 
by R = fifcL/Ls- 

Now we can rewrite the Drake equation for the observation of a CIL with microlensing and apply the numerical 
values for some the Drake factors. The result is 

Ncil,)! ^ R X (19) 

where subscript /r represent for microlensing events. For the parameter of R we take an optimistic value assuming 
that every planet that supports life will eventually develop an intelligent life (i.e. // = 1) and every intelligent life 
eventually develops an intelligent civilization (i.e. fc = !)■ For the life time of an intelligent civilization (i.e. L), while 
no signal from an intelligent civilization is detected yet, however we can assume that the present epoch of cosmology 
(after finishing the star burst epoch and present calm situation), there might be a burst of intelligent civilization in 
the Universe. So let us take L in the order of hundred million years, while there are pessimistic estimation of the order 
of 200 years (Lemarchand 2009). The hypothetical longer time scale for the life time of intelligent civilization could 
be due to developing and spreading of artificial intelligent life at larger areas of Galaxy. Then let us take L/Lg in the 
order of 10“^. Substituting numerical values for all the prameters, the number of intelligent life detection with SKA 
follow up microlensing observation would be Ncil,ii — 10“^ x iV^. 

We can estimate Ncil,^i events for two cases of single and binary leasing. Assuming 2000 microlensing observations 
per year, we would expect to have = 30 illuminating events from the binary and = 50 illuminating events 
from the single lensing where for the binary leasing Ncil,^, is in the order of 0.3 event per year and for the single 
leasing, it is in the order of 0.5 event per year. The next generation of microlensing events will certainly increase the 
number of microlensing candidates per year by probing larger areas as including the spiral arm directions and lowering 
the limiting magnitude of observation. So we may reach in near future to the number of one microlensing event with 
possible intelligent life signal from the binary lensing events (taking into account optimistic values for the parameters 
of Drake equation). 

As we noted before while the probability of detection of Intelligent civilization with the binary lensing is less than that 
of single lensing, using the early warning system, we can aware of caustic crossing of planet and as we discussed before, 
the probability of this event is about 85%. This makes the possibility of using the strategy of target of opportunity in 
the observation of intelligent life with the SKA telescope 

6. CONCLUSION 

In this work, we proposed using the gravitational microlensing as a natural telescope for magnifying signals from an 
intelligent life at the Galactic scale. Microlensing can illuminate signals of a planet by single or binary lensing mode 
for a short period of time and if in this period we observe these target with a radio telescope, these signals can be 
detected. We assumed that intelligent life is using the same wavelengths that we are using in our telecommunications. 
Telescopes as Square Kilometre Telescope (SKA) which is designed for the cosmological observation is suitable for the 
observations of these signals. 

We assumed an Earth-like military transmitter on a planet with intelligent civilization and took SKA telescope for 
detecting transmitted signals. Then we simulate the number of microlensing events with enough illuminations from the 
planets orbiting within the habitable zone that can be detected by the SKA telescope. In order to estimate the number 
of events, we performed a Monte-Carlo simulation and obtained the number of high magnification events in both single 
and binary lensing channels. We have shown that taking into account the wave optics effects in the gravitational lensing 
doesn’t change results of geometric optics. The result of simulation is that for single-lens microlensing events for 2.5% 
of overall event we will observed illumination of planet around the source star and for the binary lensing that is about 
1.5% of all the events. 

Finally we used Drake equation in our study and adapted the factors of this equation from the recent observations to 
estimate the number of intelligent life detection. For the factor of probability of intelligent life formation (i.e. //) and 
communicating civilization (i.e. fc), we adapt the optimistic values. Also we assumed that intelligent life might be 
developed in the form of artificial intelligent civilization and spread in the local areas of Galaxy. The annual number 
of events, performing follow-up observation of microlensing events with SKA telescope and getting a positive signal is 
about 0.5 for illuminating events of single lensing and 0.3 for the case of binary lensing. While the number of binary 
lensing events is less than single lensing, the probability of illuminating effect for the binary lensing is almost 85%. 
One of the advantages of binary microlensing events is that it can be identified long before the caustic crossing, using 
early warning system of surveys. On the other hand, having binary lensing events, it can be observed by SKA as target 
of opportunity mode which make this program feasible. Increasing the number of microlensing events in the direction 
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of Galactic bulge as well as spiral arms of Galaxy is achievable with the next generation of microlensing surveys, using 
wider field of view and increasing the sensitivity of detectors. This method of SETI observations would be a new item 
for the astrophysical applications of gravitational microlensing observations (Rahvar 2015). 


We would like to thank anonymous referee for his/her useful comments that helped for improving this work. 


APPENDIX 

A. WAVE OPTICS EFFECT AND SUPPRESSING OF SIGNALS 


A point like astronomical object produces a flat wave front at large distance from the source. Let us assume that 
we put a Young double slit experiment in front of this wave front. Slits at suitable separation results in fringes on the 
observer plane. The gravitational lensing system can play the same role of Young double slit experiment by producing 
double or multiple images on the lens plane . The result would be interference on the observer plane where slits are 
locates at the astronomical scales far from the observer. The wave optics property of light in gravitational lensing 
is extensively discussed in Schneider et al. (1992). This method also is a complimentary observation for detection of 
planets around the lenses (Mchrabi & Rahvar 2013). 

In geometric optics, the caustic lines on the source plane are singular lines where the Jacobian of determinant 
between the source plane and image plane is zero. However in the wave optics since the intensity of light is distributed 
in the diffraction fringes, the singularity is resolved. The magnification factor from a point like source near the caustic 
line is given by 


m(2/) 


l^max 




(Al) 


where y is the distance of source from the caustic line, Ai{x) is the Airy function and Q ~ 0.5357 is the maximum 
value of Airy function. The maximum magnification of the fringes as well as the wavelength of Airy function depends 
on the Fermat potential as 


l-^max — 


|</>ll||</222P^^ 


Yo = 



(A2) 


where (j) is the Fermat potential for a binary lens at the position of the images, subscript is the spatial derivative in 
two directions and / is the multiplication of wavenumber to the Schwarzschild radius of lens / = 2kRs- For a binary 
lenses the Fermat potential is given by 

/>(x,y) = i(x-y)2-giln(|x-xi|)-92ln(|x-X2|), (A3) 


where qi and 92 are the relative mass of lenses to the total mass of system, x is the position of image normalized to 
the overall Einstein radius, y is the position of source normalized to the projected Einstein radius of binary system 
and Xi and X 2 are the position of lenses. The lens equation from the Fermat potential obtain by extremum condition 
of this function (i.e. d(j){x,y)/dx = 0). 

Detailed calculation of Fermat potential provides |</> 222 |~^^^|</>ii|~^ — and the denominator of equation (A2) 
for equal mass lens is in the order of 0.8. For nominator of equation (A2), we adapt one meter size electromagnetic 
wavelength for intelligent life transmitter, also one solar mass for the mass of binary lenses. The result is the maximum 
magnification of planet’s signal, Umax = 75, during the caustic crossing. Using equation (8) for minimum amount 
magnification for the detection, we obtain the limiting distance of Ds < 10 kpc for observability, which is well within 
the observable distance of sources in Bulge microlensing events. Hence the wave optics formalism doesn’t change the 
results of what is calculated based on gravitational lensing in the geometric optics. 
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